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Abstract 



We investigate tlie mass spectrum of l"*" exotic mesons with quark 
5^ \ content (cscq)/(cqcs), using molecular and diquark-antidiquark op- 

erators, in quenched lattice QCD with exact chiral symmetry. For 
the molecular operator {(q7ic)(c75s) — (c7js)(q75c)} and the diquark- 
antidiquark operator {(q'^C7jc)(sC75C^) — (q'^C75c)(sC7?'c^)}, both 
detect a 1"*" resonance with mass around 4010 it 50 MeV in the limit 
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1 Introduction 

Since the discovery of Ds{2317) [1] by BABAR in April 2003, a series of 
new heavy mesons^ with open-charm and closed-charm have been observed 
by Belle, CDF, CLEO, BABAR, and BES. Among these new heavy mesons, 
the most intriguing ones are the charmonium-like states, X(3872) [3], y(3940) 
[4], F(4260) [5], Z(3930) [6], and X(3940) [7]. Evidently, one can hardly in- 
terpret all of them as orbital and/or radial excitations in the charmonium 
spectrum. Thus it is likely that some of them are exotic (non-gg) mesons (e.g., 
molecule, diquark-antidiquark, and hybrid meson). Theoretically, the central 
question is whether the spectrum of QCD possesses these resonances, with the 
correct quantum numbers, masses, and decay widths. 

Recently, we have investigated the mass spectrum of closed-charm exotic 
mesons with J^*^ = 1 [8], and 1^"^ [9], in lattice QCD with exact chiral sym- 
metry [10, 11, 12, 13, 14]. By constructing molecular and diquark-antidiquark 
operators with quark content (cqcq), we measured their time-correlation func- 
tions, and extracted the mass spectrum of the hadronic states overlapping with 
these exotic meson operators. Our results suggest that y(4260) and X(3872) 
are in the spectrum of QCD, with J^'-'' = 1 and 1"'"'' respectively, and both 
with quark content (cucu). Note that we have been working in the isospin 
limit (with rriu = frid), thus our results [8, 9] also imply the existence of exotic 
mesons with quark content (cdcd), even though we cannot determine their 
mass differences from those with (cucu). Moreover, we also observe heavier 
exotic mesons with quark contents (cscs) and (cccc), for J^'-' = l"*""*" [9], and 
1" [8]. 

Now if the spectrum of QCD does possess exotic mesons with quark con- 
tent (cqcq), then it is likely that there are also exotic mesons with other quark 
contents, e.g., (cqcs) and (cscq). However, in general, whether any combina- 
tion of two quarks and two antiquarks can emerge as a hadronic state relies 
on the nonperturbative dynamics between these four quarks. 

In this paper, we investigate the masses of the lowest-lying states (with 
J^ = 1+) of molecular and diquark-antidiquark operators with quark content 
(cscq) or (cqcs). For two lattice volumes 24^ x 48 and 20^ x 40, each of 100 
gauge configurations generated with single-plaquette action at P = 6.1, we 
compute point-to-point quark propagators for 30 quark masses in the range 
0.03 < niqa < 0.80, and measure the time-correlation functions of these ex- 
otic meson operators. The inverse lattice spacing a~^ is determined with the 
experimental input of /^r. The strange quark bare mass mga = 0.08, and the 
charm quark bare mass rrica = 0.80 are fixed such that the masses of the cor- 
responding vector mesons are in good agreement with 0(1020) and J/ijj{3097) 
respectively [15, 16]. 

"'^For recent reviews of these new heavy mesons, see, for example, Ref. [2], and references 
therein. 



2 The Molecular Operator 

In general, one can construct many different molecular operators with quark 
content (cqcs) or (cscq) such that the lowest-lying state of each operator has 
J^ = 1+. However, not every one of them has good overlap with the lowest- 
lying hadronic state having the same quark content and quantum numbers. In 
the following we only present one of the good molecular operators ^. Explicitly, 
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Figure 1: The ratio of spectral weights of the lowest-lying state of the molecular 
operator M, for 20^ x 40 and 24^ x 48 lattices at /5 = 6.1. The upper-horizontal 
line R = (24/20)"^ = 1.728, is the signature of 2-particle scattering state, while 
the lower- horizontal line i? = 1.0 is the signature of a resonance. 



We measure the time-correlation function, 



C(t)=^(M(f,t)Mt(0,0) 



^The operators we have investigated include (q7ic)(c75s), (c7is)(q75c), and 
[(q7ic)(c75s) ± (c7is)(q75c)], and they all yield compatible masses for the lowest-lying 
meson. 
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Figure 2: The mass of the lowest-lying state of M versus the quark mass TTigO, 
on the 24^ x 48 lattice at /3 = 6.1. The solid line is the linear fit. 



where the cc annihilation diagrams are neglected such that C{t) does not 
overlap with any conventional meson states. Also, C{t) is averaged over Ci 
(with 7j) for i = 1,2,3, where in each case, the "forward-propagator" Cj(t) 
and "backward-propagator" Ci(T — t) are averaged to increase the statistics. 
The same strategy is applied to all time-correlation functions in this paper. 
Then the average of C{t) over all gauge configurations is fitted to the usual 
formula 



2ma 



l^-mat _|_ g-ma(T-i)l 



to extract the mass ma of the lowest-lying state and its spectral weight 

w ' 



2ma 



Theoretically, if this state is a genuine resonance, then its mass ma and spectral 
weight W should be almost constant for any lattices with the same lattice 
spacing. On the other hand, if it is a 2-particle scattering state, then its 
mass m,a is sensitive to the lattice volume, and its spectral weight is inversely 
proportional to the spatial volume for lattices with the same lattice spacing. 



In the following, we shall use the ratio of the spectral weights on two spatial 
volumes 20'^ and 24'^ with the same lattice spacing {j3 = 6.1) to discriminate 
whether any hadronic state under investigation is a resonance or not. 

In Fig. 1, the ratio (i? = ^^20/^^24) of spectral weights of the lowest-lying 
state extracted from the time-correlation function of M on the 20^ x 40 and 
24^ X 48 lattices is plotted versus the quark mass niqa G [0.03,0.80]. (Here 
q (q) is always taken to be different from c (c) and s (s), even in the limit 
rUq —^ rric or rris). Evidently, i? ~ 1.0 for the entire range of quark masses, 
which implies that there exist J^ = 1"*" resonances, with quark contents (cscii) 
and (cscd) respectively. 

In Fig. 2, the mass of the lowest-lying state extracted from the molecular 
operator M is plotted versus rriqa. In the limit rriq -^ rriu^d — 0.00265a~^ 
(corresponding to rrij^ = 135 MeV), it gives m = 4007(34) MeV. 

3 The Diquark-Antidiquark Operator 

In general, one can construct many different diquark-antiquark operators with 
quark content (cqcs) or (cscq) such that the lowest-lying state of each operator 
has J^ = 1^. However, not every one of them has good overlap with the lowest- 
lying hadronic state having the same quark content and quantum numbers. In 
the following we only present one of the good diquark-antidiquark operators. 
Explicitly, 

D,{x) = -^ {(q'^C7,c),,(sC75C^).a - (sC7f c^),,(q^C75c),.„} (2) 

where C is the charge conjugation operator satisfying C'y^C~^ = — 7J and 
(6*75)^ = — C75. Here the "diquark" operator {q^TCD^a for any Dirac matrix 
F is defined as 

(q TQjxa = €abcq.xab^al3Q,xl3c (3) 

where x, {a,b,c} and {a,P} denote the lattice site, color, and Dirac indices 
respectively, and eabc is the completely antisymmetric tensor. Thus the diquark 
(3) transforms like color anti-triplet. For F = C75, it transforms like J^ = 0^, 
while for F = C^i {i = 1, 2, 3), it transforms like l"*". 

In Fig. 3, the ratio {R = W20/W24) of spectral weights of the lowest-lying 
state extracted from the time-correlation function of D4 on the 20^ x 40 and 
24^ X 48 lattices is plotted versus the quark mass rriqa G [0.03, 0.8]. Evidently, 
i? ~ 1.0 for the entire range of quark masses, which implies that there exist 
J^ = 1+ resonances, with quark contents (cdcs) and (cues) respectively. 

In Fig. 4, the mass of the lowest-lying state of the diquark-antidiquark 
operator D4 is plotted versus niqa. In the limit m^ -^ rnu^d, it gives m = 
4015(25) MeV. 
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Figure 3: The ratio of spectral weights of the lowest-lying state of diquark- 
antidiquark operator D^, for 20^ x 40 and 24^ x 48 lattices a.t (3 = 6.1. The 
upper-horizontal line R = (24/20)^ = 1.728, is the signature of 2-particle 
scattering state, while the lower- horizontal line i? = 1.0 is the signature of a 
resonance. 



4 Summary and Discussions 

In this paper, we have investigated the mass spectra of the molecular operator 
M and the diquark-antidiquark operator D4, in quenched lattice QCD with 
exact chiral symmetry. Our results for M and D4 are summarized in Table 
1, where in each case, the first error is statistical, and the second one is our 
crude estimate of combined systematic uncertainty. 

Evidently, both the molecular operator M and the diquark-antidiquark 
operator D4 detect a l"*" resonance around 4010 ± 50 MeV in the limit niq — >■ 
rriu^d- Since its mass is just slightly above Y(3940), high energy experiments 
should be able to see whether such a resonance, say, Xg, exists in some decay 
channels, e.g., Xg -^ KttJ/'^, in the near future. 

Note that in our calculations, we have omitted all internal quark loops 
generated by vacuum fluctuations (i.e., the quenched approximation). How- 
ever, we expect that the quenched approximation only results a few percent 
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Figure 4: The mass of the lowest-lying state of the diquark-antidiquark oper- 
ator D4 versus the quark mass rriga, on the 24^ x 48 lattice at /3 = 6.1. The 
solid line is the linear fit. 



systematic error in the mass spectrum, especially for charmed mesons [16]. If 
it turns out that high energy experiments do not find any charmonium-like 
resonances around 3960-4060 MeV, in particular, in the KttJ/'^ channel, then 
it might imply that quenched QCD could be dramatically different from the 
unquenched QCD, even for the mass spectrum. 

Acknowledgement 

We thank Steve Olsen for his question regarding (cues), which motivated 
our present study. This work was supported in part by the National Sci- 
ence Council, Republic of China, under the Grant No. NSC94-2112-M002-016 
(T.W.C.), and Grant No. NSC94-2119-M239-001 (T.H.H.), and by the Na- 
tional Center for High Performance Computation at Hsinchu, and the Com- 
puter Center at National Taiwan University. 



Operator 


Mass (MeV) 


R/S 


^2 (U7ic)(c75s) (C7is)(u75c)] 
^ {(u^C7.c)(sC75C^) - (u^C75c)(sC7fc^)} 


4007(34) (31) 
4015(25)(27) 


R 
R 



Table 1: Masses of the lowest-lying states {J^ = 1"*") of the molecular operator 
M and the di quark- antidiquark operator D4. The last column R/S denotes 
resonance (R) or scattering (S) state. 
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